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NATIONAL AD7IS0RY CCMMITTSE FOR AERONAUTICS 


TECHNICAL NOTE NO. 1343 


CRITICAL STRESS OF THIN-V'ALLED CYLINDERS IN AXIAL COMRRESSION 
By S. B. Batdorf, Murry Schildcrout, and Manuel Stein 


SUMMARY 


Empirical design curves are presented for the critical 
stress of thin- walled cylinders loaded in axial compression. 
These curves are plotted in terms of. the nondlmenslonal 
parameters of small-deflection theory and are compared with 
theoretical curves derived for the Buckling of cylinders with 
simply supported and clamped edges. An empirical equation is 
given for the Buckling of cylinders having a length-radius 
ratio greater than aBout 0 .75 • 

The test data oBtained from various sources follow the 
general trend of the theoretical curve for cylinders with 
clamped edges, agreeing closely with the theory in the case of 
short cylinders, But falling consideraBly Below the theoretical 
results for long cylinders. The discrepancy in the case of 
long cylinders increases with increasing values of the ratio 
of radius to wall thickness. Plotting curves for different 
values of this ratio reduces the scatter in the test da,ta 
and a certain degree of correlation with theory is achieved. 
Advantage is taken of this correlation to oBtain estimated 
design curves for cylinders with simply supported edges, for 
which little experimental information is availaBle. 


REOT'T OF PRS7I0US WORK ON PROBLEM 


Solutions to the proBlem of the determination of the critical 
stress of thin-walled cylinders suBJected to axial compression 
have Been presented By a large numBer of authors. Southwell, 
Timoshenko, Flugge, and numerous other authors have oBtained 
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theoretical solutions "by the use of the small-deflection theory. 
(See, for example, references 1 to l^-.) The value given "by the 
small-deflection theory for the 'buckling stress of a thin-walled 
cylinder of moderate .length having simply supported edges is 


0.608 a 

Wid - m5) , 


( 1 ) 


fdiere 

C3:*'''tical compressive stress 
E Young* s moduliis ' . - 

t wall ttiickness of cylinder 
r radius of cylinder 

H Poisson's ratio (in the pi-esent paper p is taken to be 0,316 
whenever a value is assigned to it) 

Experiments (references 5 to 10) have shom that the ^tu^ 
critical stress is much lower than that predicted "by equation (l). 
Except In the case of short cylinders, the experiments usually 
give values only I 5 to 50 percent cf that predicted theoretically; 
moreover, the observed buckle pattern is different from that 
predicted on the basis of theory. A number of attempts have been 
made to explain these discrepancies theoretically, Elligge (refer- 
ence 3 ) considered tlie deviation of the actual edge supports , from' 
the support conditions assumed in the theoretical treat^nt. Donnell 
(reference 5 ) and also Eliiggo considered the initial deviation from 
the perfect cylindrical shape. Neither of the two attempted 
exp.lanatlons satisfactorily accounts for the d.lscrepancy existing 
between the theoretical .and expei’imental values of .the buckling 
stresses of cylinders. 

"V'on Kajnaan and Tsien (reference .11) introdiiced a lai’ge- 
def lection 'theory to acco\?nt for the buckling behavior of long 
cylinders. They showt3d that a long cylinder can be in equilibrium 
in a bxickled state at a stress , that is much smallo-r than the 'critical 
stress of small-deflection theory and also succeeded’ in accounting 
for the buckle pattern observed in the early stages of buckling. 
Reference 11 suggested that idien a cylinder has an initial imperfec- 
tion or is subjected to a shock, it might pass into one of these 
buckled states without ever having reached the critical load given • 
by equation (1). Based on the same approach, a theory for tho 
buckling stresses of perfect cylinders was proposed by Tslen 
(reference 12), which gave 
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°x ='0’370 ~ (loading liy rigid screw -power testing machine) (2) 


= 0.238 (loading hy , ideal hydraulic testing machine or 

^ dead weight) (3) 


The large -deflection theories fall in two respects to descrihe 
fide^uately the buckling behavior of actual cylinders. First, the 
theories are formulated only for long cylinders; eguations (2) and (3) 
seriously underestimate the critical stress of very short cylinders. 
Second, even for long cylinders, attempts to determine eaperi- 
mentally the numerical coefficient C in the buckling formula 


a Cfii 

X 


w 


have resulted in appreciable experimental scatter. The experimental 
scatter is due at least in part to the initial imperfections of 
construction always present in real cylinders. (See fig. 1 taken 
from reference I3, fig. 18.) 

In the, absence of a complete and satisfactory theoretical 
solution for the critical stress of cylinders, a number of authors 
have proposed empirical foimilae derived from test data (refer- 
ences 6 to 8) ., One such formula, which takes into accomt the 
length of the cylinder, is due to Ballerstedt and Wagner (reference 8) 


X = 3.3(i)" 0.2 i (57 

The first parameter in this equation is appropriate for flat 

j, \L/ 

sheet and the second paraneter i is included to take into account 

the effect of curvature. More recently Kanemitsu and Nojlma 
(reference 9) ccmpiled all available previous experimental results 
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and conducted a number of teste of their o\m. The formula of 
Weigner and Bailer stedt was modified in reference 9 "to bring it into 
better agreement with experiment as follows: 


'X 


E 




Within its range of application 




O.K ^<1.5; 500<f< 3000 


■) 


(€) 


equation ( 6 ) is in considerably better agreement idth experiment 
than equation ( 5 ) but, because of the change in the exponents of 
the parameters, eqtiatlon ( 6 ) does not have any rational basis and 
must be regarded as purely empirical. A completo divorce of theory 
and experiment, however, cannot be regarded as a satisfactory 
permanent settlement of the problem, and the present report attempts 
to bring theory and experiment into jroasonable accord. 


CONTRIBUTION OF ERESENT PAPER 


In the present paper the available test data for critical 
stresses of cylinders are reexamined, and theoretical results are 
used as a guide in fairing the curves, in extending the range of 
validity of the existing empirical results, and In achieving a 
more rational Interpretation of the test data. For this pmpose 
the test data are plotted in terms of the parameters of cylinder 
theory and are compared with theoretical results derived in the 
appendix on the basis of small-deflection theory. 

The cylinder- theory parameters used are 
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where 
D 


L 

Z 


Et^ '' 

flexural stiffness of plate per unit length 1 — 

\ 12(1 - ix% 

length of cylinder 


cuiVature parameter 

critical stress coefficient appearing in the equation 


cr = 


L^t 


The e^qierimental data are used as the principal guide in 
determining the critical campressire stresses of long cylinders 
(large values of Z) and the theoretical results are used mainly 
to supplement the test data in determining the critical stress 
of very short cylinders ( small values of Z) . The experimental 
scatter is reduced hy. presenting different curves for cylinders ; 
with different values of the ratio of radius to wall thickness 
on the asBumption that for' long cylinders this ratio furnishes some 
indication of the initial imperfections of the cylinder. Althou^ 
these curves were determined partly on the hasls of theoretical 
considerations, they are for convenience referred to herein as 
empirical curves. 


. EESULTS AND CONCLUSIONS 


The critical compressive stress for -cylinders is given hy • 
the eq.uation 


L^t . 


( 7 ) 


where the values of ^ obtained from figure 2 for 

cylinders with either clamped or simply supported edges. TVip> 
design curves for cyliijders with clamped edges are es-bahlished 
hy -the test results .reported in references 5 to 9. (See fig. 3.) 
Each curve v/as faired throu^ a series of test points which were 
plotted for cylinders with nearly the same ratio of radius to wall 
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thickness r/t. The estimated (dashed) parts of the design curves 
for simple support were obtained by fairing between the known 
experimental curves for long cylinders (large values of Z), which 
according to theory should be the same whether the cylinders have 
simply supported or clamped edges, and the theoretical curves for 
very short cylinders (small values of Z) . 

For long cylinders the buckling stress is considerably below 
the theoretical buckling stressj the amount of the discrepancy 
depending on the ratio of radius to wall thickness. For very short 
cylinders the values of the critical stresses approach those for 
flat plates (simply supported ends, k^^ = 1; clamped ends, 
kj. = 4), for which the agreement between theoretical and experimental 
results is known to be good. The general trend of each empirical 
curve is similar to that of the theoretical curve, indicating the 
existence of a certain degree of correlation between theory and 
test data. 

At large values of Z the curves for k^ become straight 
lines given by the formula 

• kj-^l.l^CZ (8) 

where C depends on the ratio of radius to wall thickness of the 
cylinders in the manner shown in figure 4. From equations (?) 
and (8) the following expression for the critical stress is obtaiiied 


■ CTj * ' ^9) 

Equations ( 8 ). and ( 9 ) may be used when the length of the cylinder 
is more than about 3/4 of the radius. The omplrical curves of 
reference 10 indicate that the critical stress is substantially 
independent of length when the length is greater than about 3/4 of 

the radius. fThis result may be checked by noting that for Z > 0.5 

the experimental curves of figure 2 are substantially strai^t lines 
of unit slope . ) 

Tn figure 5, the empirical formula of Kanemltsu and Nojima 
(equation (6) of the present paper, and the best previously published 
formula for the biAckling of cylinders) is plotted in terms of the 
parameters k^ and Z. The curves are cut off at those values of Z 
corresponding to the lower limits of the range of dimensions within 
which the formula was Intended to apply. In genereQ._, for the range 
covered, the curves are in .reasonable- agreement with the test data 
and with the curves of the present paper for cylinders with clamped 


cH*-S 
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edges. The practical importance of the present approach lies in the 
fact that the use .of - the theoretical parameters and the theoretical 
solutions as a guide in fairing the curves permits the resnoval of 
the. lower limits on these curves and also permits estimated curves 
to he dram for the huckling stresses for simply supported cylinders 
althou^ esi)erimental data are availahle only for cylinders with 
clamped edges. 


Langley Memorial Aeronautical Laboratory 

National Adviso:ry Coiamittee for Aeronautics 
Langley Field, Ya. , March 20, 191^7 . 
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APPENDIX 

THEORY POE CYLINDERS BUCKLUTG UNDER 
AXIAL COMPRESSION ■ 
Symbols 


m positive integer 

r radius of cylinder 

t wall thlclsness of cylinder 

V radial component of displacement, positive outward 

X axial coordinate of cylinder 

y circumferential coordinate of cylinder 

"fc 

C coefficient appearing in = CE - 
D flexural stiffness of plate per unit length 

E Young's modulus 

L length of cylinder 

Q operator on v defined in appendix 



Et' 


12 (1 



«m 

^x 




coefficient of deflection function 

I 

critical-ccairpresBlve- stress coefficient appearing in -ttie 

g 

formula c„ » ^ 

^ L^t 

7?To - fm - 1)^1. 

X 


r; -2p i2z^(m - i>^ 

|(m - 1) H. sg . '2 


deflection function defined in the appendix 
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P ^ 
X 


X half wave length of hiicld.es in circumferential direction 
[1 Poisson* a ratio 

critical axial compressive stress 

a" . 


^ a . 


2 a^ +.a^, 


ax^ ax^^ay^ ap* 

V”^ the inverse of defined hy V^V^v = w 


Theoretical Solution 


The critical compressive stress at which hucld-ing occurs in 
a cylindrical shell may he obtained hy solving the equation of 
equllihrlum. 


Equation of equilihrlum. - The equation of equllihrium for 
a aUghtiy huft yift d cyli ndr ical shell under axial compression is 
(reference i4) 


D 


V^' 


w + V 


rh a^w 


+ cr^t 


a^w 

as^ 


(Al) 


where x la the coordinate in the axial direction- and y is the 
coordinate in the circumferential direction. The accompanying 
figure shows the coordinate system used in the analysis. 
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Diriding throng egiiation (Al) "by D gives 

12Z^ hK . ^2^ 


V, w + 


■where the dimensionless parameters Z and Ic^ are defined by 


(A2) 


z = 


rt 


k = 

^ Dn2 


The eqmtion of eq.villit»rim may be represented by 

Qw s= 0 


■where Q is defined by 


Q = V 


k 


+ 


12zf -h 



(A3) 


Method of solution .- Eq.\iation (A2) may be sol^ved by use of 
■the Galerkln rae-fchod as outlined in reference 15. When ■this method 
is applied, ■fche deflection v is esi>ressed in series form as 
follows 


w = 



(Ah) 


The set of functions are chosen to satisfy the boundary 
conditions but need not satisfy -fche eqmtlon of equilibrium. The 
coefficients ^ are determined by the equations 

2K L ' ' ' 

/ P V- QwdxxJy = 0 (m = 1,2,3, . . . j) (A5) 

‘'0 Tt 

In •fche, present paper the deflection functions .■were chosen •fco 
satisfy the follo^wlng conditions on v at the ends of the cylinder: 
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For simply supported edges 


V t= 



= 0 


For clamped edges 


w= 0 


Simply- supported edges.- An. espansion. for vr that is sinusoidal 
in -Uae ' circimif erenila^ dli?eotion and perfectly general (subject to 
the boundary conditions for simple support) in the axial direction is 


V = sin 



5~ sm sin 


mroc 

L 


(a6) 


•vdiere X is the half ifave length of the buckles in the circum- 
ferential direction. Egoiation (a 5) is egiiivalent to equation (a 4) if 


= sin ^ sin — 
XL 


(at) 


Substitution of expressions (a£) and (AT) into equation (A5) and 
integration over the limits indicated give 


m^kx 



+ m‘ 


,2^2 


12Z^^ 
Jt^(m^ + p2^2 


(A8) 


■where 


• P = ^ = 1. 2, 3,--.) 

The minimum value of k^^ for a given S is found by assuming 
a value for m and minimizing kj- -wi-th respect -to p. This pro- 
cedure is followed for various values of m until a minimum k^ is 
reached. Figure 2 presents the -theoretical critical stress coeffi- 
cients for cylinders -with siiiply supported edges subjected to avT 
compression. 


Clamped edges .- A procedizre similar "bo -that used for cylinders 
with simply supported edges may be follo-wed for cylinders -with 
clamped edges. The deflection function used is the folio-wing series 


V = sin 


is: v“ 

X 2— 


/ 

aia^os 


(m - l)rtx 


- cos 


(m + l)rtx 




(A9) 
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Each term of ■tiois series satisfies the conditions on v at the 
edgss. The fijncticn is nov defined as 




( m - l)]Qc 

L 


cos 


(m + l)Tfx 




(m « 1, 2,. 3, . . .) (AlO) 


After the same operations are carried out for clamped edges as 
those carried out for the caae of simply supported edges, the following 
eouatidns result: 

For m *= 1 

ai (2Mi + M 3 ) - SgMg = 0 

For m = 2 


Sg (Mg + M 4) — a^M^ ® 0 


For m = 3, 4, 5, . . . ' 


> 


(All) 


®m %+s) ~ " ®m+2^+2 = 0 


where 




&-1) 


2 ^ p2 


12Z^(m - 1)^ , .0 

. (a - 1)2 


j{ni - 1)^ + (30 


(m = 1, 2, 3> • • • ) 


These equations have a solution if the following infinite determinant 
vanishes: 
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m=l 

2 Mi+I 43 

0 

-M 3 

0- 

rte2 

0 

Mg-rtlit 

0 

-M4 

nts3 

-M3 

0 

M^4M- 
-> J 

0. 

Uteit 

0 

-M4 

0 


2 h =5 

0 

O' 

-M 5 

0 

m=6 

0 

0 

0 


nt=7 

0 

0 

0 

• 0 

EfcS 

0 

0 

0 

0 


0 

0 

0 

0 . 

0 

0 

0 

0 


0 

0 

0 

0 

-% 

0 

0 


0 


0 

0 

Mg 4 Ma 

0 


i 

0 

M^ 4 M^ 

0 

0 

-1% 

0 
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(A 12 ) 


If the rows and coliartns are rearranged the infinite determinant can 
he factored into the product of two infinite siihdetoi’minants . The 


cri- 

bicaJ. stresses may 

then he 

obtained from the foUovlng eq.u 

atlon 

13=1 

2Mi-iM3 

-M3 

0 

0.. . 

0. 

P 

0 

0 ... 


tfe =3 

-M3 

M3-tM^ 

-M^ 

0... 

0 

0 

0 

0 ... 


ei =>5 

0 

-M5 


—My • • « 

0 

0 

0 

0 ... 


ib =7 

0 

0 

-‘V 


0 

0 

0 

0 ... 


« 

• 

« 


• 

• 

• 

• 

» • 



f 

• 

* 

• 

• 

* 

» 

m 


• 

• 

• 

« 

• 

9 

9 

• 

• 

= 0 

ri=2 

0 

0 

0 

0... 



0 

0 ... 


Efc =4 

0 

0 

0 

0.. . 

-M^ 


-% 

0 ... 


za =6 

0 . 

0 

0 

0. . . 

0 



-Mg ... 


m=8 

d 

0 

0 

0. . 

0 

0 

• 

-Mg 

%"^1Q • • • 


« 

• 

• 

• 

• 

• 

• 

• 

• 


• 


« 

• 

• 

• 

9 

» 

* 


• 

• 

• 

■ 


9 

9 

• 

9 

• 



The infinite euhdetermlnant involring terms with odd suhscripts 
corre spends to a symmetrical huclding’ pattern (a hucKLlng' pattern 
synmetrical about the plane perpendicular to and hieecting the axis 
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of the cylinder). The infinite suMeterminant Involvlixg terms vlth 
even subscripts corresponds to an anti symmetrical buoKLing pattern. 
For brevity these subdeteamilnants •will be referred" to as the odd 
determinant and the even determinant, respectively. 


The first approximation is 


Odd determinant; 

2 M 2 -+ = 0 

(Al4) 

Even determinant; 

M 2 + M]^ e 0 

(A 15 ) 


The second approximation is 

Odd determinant; + m^) + ^ 3 ^ = 0 . ■ 

Even determinant; M 2 (Mi^ + Mg) + MijMg =* 0 (AI 7 ) 


These eq.uations show that for a selected value of the curvature 
parameter Z the critical buckling stress of a cylinder depends 
upon the circumferential t^ave length. Since a struct'ure buckles at 
the lowest stress at -which instability can occur, kx is minimized 
■with respect to the -wave length by substi'fcuting valties of 3 into 
the eq-uations until the minimum can be obtained from a plot 
of kx against 3. For a given Z the lower of the t-wo values 
obtained from equations (Al4) and (AI5) is the first approximation 
of the critical buckling stress and, similarly, the lower of the 
ttro values of kx obtained from equations (A16) and (AI 7 ) is the 
second approximation of the critical buckling s-tress. 

Figure 2 presents the theoretical critical stress coefficients 
for cylinders •vd.th clamped edges in axial ccmpression as obtained 
from the second approximation, -together -wi-th -the exact curve for 
the case of simply supported edges. Althou^ this solution is an 
upper-limit solution, the second approximation for the critical 
stress coefficient of a cylinder wi-th clamped edges must be very 
close -to being exact for intermediate and large values of Z 
because it is almost iden-uical with the exact solution for a 
cylinder wi-th simply supported edges, and -the critical stress of 
a cylinder with clamped edges cannot be less than the critical 
stress for a cylinder with simply supported edges. For values of 
Z approaching zero, the accuracy of the second approximation is 
Indicated by -the fact that it coincides with the known exact 
solution (kjj. a 4) for a long flat plate with clamped edges. 



NACA OT No . I 3 U 3 


15 


REFERENCES 


1. Southwell, R. V.; On the General Theory of Elastic Stability. 

Phil. Trans. Roy. Soc. (London), ser. A, vol. 213 , 191i^, 
pp. 18 t-241<-. 

2. Timoshenko, S.: Theory bf Elastic Stability. MqOraw-Hill Book Co., 

Inc., 1936 , ch. IX- 

3- PlUgge, W.: Statik und Dynamik dor Schalen. Julius Springer 

(Berlin), 193^^, pp- 1^-199* 

k. Dean, V. R.: On the Theory of Elastic Stability. Proc. Roy. 

Soc. (London), ser. A, vol. 107, no. April 1, I 925 , 
pp. 73 ^^- 760 . 

5* Donnell, L. H.; A New Theory for the Buckling of Thin Cylinders 
under Axial Compression and Bending. A.S.M.E. Trans., 
vol. 56 , no. 11, Nov. 1934 , pp. 795 - 806 . 

6 . Lundqulst, Eugene E.; Strength Tests of Thln-Walled Duralumin 
Cylinders in Compression. NACARep. No. 473, 1933 • 

7« Wilson, Wilbur M., and Newmark, Nathan M.: The Strength of 

Thin Cylindrical Shells as Columns. Eng. Exp. Sta. 

Bull. No. 255 , Dhiv. 111., Feb. 1933 . 

8 . Ballerstedt, W., end Wagner, H.; Yersuche .uber die Festlgkeit 
dhnner unversteifter Zylinder unter Schub- und Langskr^ten 
Luftfahrtforschung, Bd. I 3 , Nr. 9 , Sept. 20, 1936, 

pp. 309 - 312 . 

9* Ean633iltsu, Sunao, and Nojima, Noble M.J Axial Compression Test 
of Thin Circular Cylinders. Calif. Inst. Tech. Thesis, 1939- 

10. Robertson, Andrew; The Strength of Tubular Struts. R. &. m. 

No. 1185 , British A.R.C., I 929 . 

11. von Earman, Theodore, and Tsien, Hsue-Shen: The Buckling of 

Thin Cylindrical Shells under Axial Compression. Jour. 

Aero. Sci., vol. 8 , no. 8 , June 1941, pp. 303 - 312 . 

12. Tsien, Hsue-Shen: A Theory for the Buckling bf Thin Shells. 

Jour. Aero. Sol., vol. 9> no. 10, Aug. 1942, pp. 373-384. 



16 


MCA TN No. 1343 


13 . Donnell, L. H<>; The Stability of Isotropic or Orthotropic 

Cylinders or Flat or Curved Panels, between and across 
Stiffeners, with Any Edge Conditions between Hinged and 
Fixed, under Any Combination of Ccmpression and Shear. 
MCA IS No. 918 , 1943 . ■ 

14. Batdorf, S. B. : A Simplified Method of Elastic-Stability 

Analysis for Thin Cylindrical Shells. H - Modified 
Equilibrium Equation. NACA IN No. 1342, 1947. 

15 . Duncan, W. J,; The Principles of the Galerkin Method. 

E. & M. No. 1848,. British A J?.C., 1938. 



"Perfect" strut 




NATIONAL ADVISORY 
COMMITTEE RM AEMMAUTICS 


Figure 1,“ Effect of Initial defects upon the maximum 
compressive loads of struts and cylinders. 

(From reference 13, fig, 18.) 
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Figure 3.- Comparison of test re 
recommended for cylinders wit 
(Data from references 5 to 9. 
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